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Abstract 

In the standard model, energy release in dense stars is severely restricted by the 
Pauli exclusion principle. However, if, in regions of space of high fermion degeneracy, 
there is a phase transition to a state of exact supersymmetry (SUSY), fermion to 
sfermion pair conversion followed by radiative transitions to the Bose ground state 
could lead to a highly collimated gamma ray burst. We calculate the cross section 
for ee — ► ee in a SUSY bubble and construct a monte carlo for the resulting sfermion 
amplification by stimulated emission. 



1 Introduction 

Recently it has been proposed that the puzzling gamma ray bursts seen over the last few 
decades in satellite observations pQ could be due to a transition to the SUSY vacuum in 
compact astrophysical objects such as white dwarf or neutron stars [2]. This phase transition 
could be catalysed by the high matter density of such stars and, once nucleated, will spread 
through the entire star. Such a catalysis has been demonstrated in low dimensional models 
[3] at high density and is probably a feature of dense matter in arbitrary dimensions [3]. 

In such a SUSY bubble, normal particle pairs can convert to sparticle pairs, which, being 
bosons, can drop to the ground state through gamma ray emission. In this way, the entire 
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kinetic energy of a degenerate fermi sea can be radiated away. It is possible that the high 
degree of collimation observed in the gamma ray bursts could be due to Bose enhancement 
as is familiar in terrestrial lasers. If this transition occurs in white dwarf stars, there will 
be a multi-component structure to the bursts as various particle species convert to their 
SUSY partners, possibly interrupted by periods of fusion of super symmetric nuclei. A more 
complete discussion of the physical picture proposed is contained in [2]. This model is 
dependent on the assumption that the common particle and sparticle mass in the exact 
SUSY phase is equal to (or less than) the particle mass in the broken SUSY phase. This 
mild though necessary assumption is, perhaps, supported by the string theory result that 
the ground state super mult iplets are of low (in fact zero) mass. 

The possibility of a phase transition between vacua requires that the vacuum structure 
of the theory is dynamically determined as in string theory and in certain other models of 
susy breaking. The other possibility, that the susy breaking is determined by arbitrary fixed 
parameters as in the minimal supersymmetric standard model is probably less satisfying from 
a theoretical perspective. The transition we consider (from an unstable de Sitter vacuum 
with positive vacuum energy to a stable susy vacuum) has been considered in the formal 
string theory study of ref.jH] although phenomenological consequences are not part of that 
work. 

In this article we treat the simplest component of the phase transition, namely the con- 
version of electron to selectron pairs. 



This cross section was calculated previously in |6| neglecting the electron mass as appropriate 
in the broken SUSY phase where the electron mass is many orders of magnitude less than 
the selectron mass. In the exact SUSY phase the electron and selectron masses are equal 
and the cross section near threshold is needed. The corresponding simple extensions of the 
cross section formulae are given in section II. The results are applicable to a possible SUSY 
conversion in a white dwarf which, in the broken SUSY phase, is stable against gravitational 
collapse due to the electron degeneracy. In a neutron star, the SUSY conversion would 
be more complicated and the energy release would be more slow since sneutrons will not 
efficiently radiate. Even in the white dwarf case, nuclear processes may be somewhat more 
important than the electron component which we treat here. In section III, we use the cross 
section of ll.ll in an event generator for the differential transition rate incorporating the Bose 
enhancement effect. We assume that the entire kinetic energy of the selectrons is released 
into gamma radiation as the scalars drop into the ground state although we do not treat 
the radiative processes explicitly in this paper. Section IV is reserved for a summary of our 
conclusions. 





2 Cross section calculation 



In the process e e 
in figure 1. 



ee the mediator of the interaction is a Majorana spinor 7 as shown 



We follow the Feynman rules as reviewed in |7]. In particular the vertex factor for an 
incoming left (right) handed Dirac fermion of spinor index /3, an outgoing SUSY partner, 
and an outgoing Majorana fermion of spinor index a is 
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(2.i; 



and for a Majorana fermion line with double incoming arrows as in figure 1 carrying mo- 
mentum q from this vertex to one where the Majorana fermion is emitted with index 7, the 
propagator is 
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q 2 — m 2 + ie 

Acting on the Dirac spinor the charge conjugation operator, C, has the effect 

u T (p, s)C~ 1 = —v(p, s). 



(2.2) 
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For the associated production of left and right selectrons, the amplitudes corresponding 
to graphs a and b in figure 1 are then 
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Here and throughout the paper the Mandelstam variables are 



s=( Pl +p 2 ) 2 
t=( Pl -p 3 ) 2 
u=(pi - Pa) 2 . 

The squared matrix elements take the form: 

\M LR \ 2 = £ ss , (M a Ml + M h M\ + M a M\ + M h M\ 

ut — 2tm 2 e — (jnl — rr 

Am 2 (ml — in 2 
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where S ss / denotes averaging over the spins of incoming particles. The production of two 
left or two right selectrons is obtained from the above by the appropriate changes in the 




Figure 1: Feynman graphs for the conversion of an electron pair to selectrons, eq 11.11 via 
photino exchange . 

helicity projection operators and dividing by the statistical factor for the identical final state 
particles. 

\M RR \ 2 = \M LL \* = ^(s- 2ml) { t\t, + -A^) 2 • (2-5) 



2! V~ --ey \t-M? u-M? y 

In the limit of negligible electron mass, these agree with the spin-averaged formulae of [7j 
and with the squared helicity ampltudes of [Hj. Other authors have considered the effect of 
other neutralino exchanges jSJ El but, since these particles in the SUSY phase have masses 
comparable to the W and Z, they are negligible for our purposes. 

In the numerical calculations of Section III, we will need the integrated matrix elements 
squared: 

fA B (s)= i 1 /„ dt\M AB \ 2 . (2.6) 

Js(s — Am 2 ) 
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The total cross sections are related to the Jab by 



a T {ee -> e A e B ) = (2.7) 
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Calculations for the different Matrix elements in the SUSY phase (m e = m g = m) will 
give us the following: 
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Figure 2: Behavior of the t integrated matrix elements squared as a function of s — 4m 2 . 

These are only logarithmically sensitive to the photino mass. The selectron momentum 
and angular distributions to be calculated in the next section are even less sensitive to the 
photino mass since they depend only on the shape of the cross sections and not on the 
absolute values. Integrals of Jll and Jlr as required in the conversion rates discussed in the 
next section are finite in the limit of zero photino mass. However, a more precise treatment 
of photino mass effects requires a careful treatment of experimental resolution beyond the 
scope of this study. The behavior of the functions /ll(s) and }'lr{s) for My = m/4 is shown 
in figure 2. For larger values of the photino cutoff mass, f LR lies below f LL although the 
near constant values of the two f's well above threshold are not affected. 



3 Event generation 



In a typical white dwarf (solar mass and earth radius) there are 

iVo = 6 • 10 56 



(3.1) 



electrons in a degenerate Fermi sea: 



M = **. (3.2) 



The Fermi momentum is 



Pr, 



3N \ 1/3 



2nh = 0.498(— ) 1/3 MeV/c. (3.3) 
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Here n is the electron number density and n is that of the nominal white dwarf of solar 
mass and earth radius. The event rate for process ITTI in a volume V is given in terms of the 
differential cross section by 

T = f d a{^yN l . (3.4) 

where the incident and target distributions are given bv 13.21 The differential cross section is 
given in terms of the invariant matrix element squared by 

doJ^Jf 1 '. (3.5) 
AE x E 2 v v ; 

Here m is the common electron and selectron mass in the SUSY bubble. The energies, E{ 
and the incident velocity, v, are those of the target rest frame. The final state phase space is 

= fP l d^ ± 5\p 1 +p 2 -p 3 - P4 ) 
s 2E 3 2E 4 (2tt) 2 v ' 

where any statistical factors that occur for identical final state particles are put into the 
matrix element squared. Thus the Lorentz invariant event rate per unit volume is 

dY 4| M | 2 d 3 pi d 3 p 2 d 3 p 3 d 3 p 4 4 

V = ^^^2E 2 2E 3 2E A 5 iPl + P2 ~ P3 ~ Pi) - (3 ' 7) 

The P2 integral can be done trivially using the 5 function and the pi integral is then conve- 
niently done in the CM frame with p 3 = e z . 

d 3 pi sf/ 2 2 ndt 
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Thus 
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or 
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For p max we use the value of E31 corresponding to the Fermi energy in a white dwarf of earth 
radius and solar mass. 

The final state of the process consists of two distinct species of scalars, Bl and Br. Thus 
the effective matrix element squared in 13.101 and elsewhere above is actually 

I M | 2 = | M LR | 2 + (| M LL | 2 + | Mrr I 2 ). (3.11) 

In a bath of pre-existing selectrons each matrix element squared is related to the elementary 
matrix element squared calculated with no pre-existing selectrons by the Bose statistical 
factors 

I M | 2 =| M 0LR I 2 {N L {p 3 ) + l){N R (p\) + 1) + (| M 0LL I 2 (N L (p 3 ) + l){N L {p,) + 1) 

+ | M 0RR I 2 (N R (p 3 ) + l)(iV B (p 4 ) + 1)). (3.12) 

The matrix elements of the previous section, calculated for the case of no selectrons in the 
initial state, are the A4q of this section. For the present we treat only the LL final state. 
The RR final state implies that there could be at least four gamma ray jets in each burst and 
possibly many more if different energy levels in the Fermi sea lead to independent gamma 
ray bursts. This suggests a picture where the gamma ray jets are much more numerous in 
each stellar explosion and, individually, much more narrow and less energetic than currently 
assumed. A more detailed analysis, beyond the scope of the current paper, is needed to 
explore this point but we speculate that this could be the cause of the "spikey" nature or 
rapid time variability of many of the observed bursts. 

Our remaining analysis in this paper, restricted to the LL jets, allows us to drop, the L 
subscript on the occupation numbers. The CM energy a/s is determined by p% and p±. The 
t integral has been done analytically in section 2: 

We define a grid of Nf, in points in each of the six variables. In practice we choose 
N bin = 15. The discretized i'th integration variable in l3.1Ul is 

Vi = Vi,min 4~ [pi,max \N% 4" l/2)/A^j n (3.13) 

where 

< Ni < N Un - 1. (3.14) 

In order to generate events with probability defined by 13.101 and 13.121 it is convenient to 
linearize the six dimensional integral. We define the composite integer variable 

8=1 

with limits 

0<j<iVL-l. (3.16) 

Each value of j corresponds to a unique value of each of the six integration variables. The 
distribution of events is then defined as an integer valued array with N^ in grid points. To 



(3.15) 



handle 15 6 grid points requires careful memory management techniques. * The integer j 
can be decomposed into two integers j 3 and j\ which encode the three dimensional vectors 
•pz and pi respectively. 

j 4 =j mod Nl n (3.17) 

The first of these equations is defined by integer division, i.e. js is the largest integer less 
than or equal to j/N^ in . The event generation follows standard techniques [TQ] except that 
the probability distribution changes with each event due to the Bose enhancement factors. 
In the event generation of the selectron distributions, constant overall factors in I3.1UI play 
no role. We begin by putting all the selectron occupation numbers to zero and calculating 
the (unnormalized) probabilities 

P(j) = hL(s)^r(N(j 3 ) + l)(N(j 4 ) + 1). (3.18) 



We then define the partial sum 



i=0 



as well as the complete sum 
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E m (3-20) 



i=0 



R(j) is a monotonically increasing function of j. One then calls a random number w between 
zero and one. For the unique value of j for which w > R(j — 1)/Pj nt and w < R(j)/Pi nt one 
increments the selectron occupation numbers by one: 

N( j3 ) - N{j 3 ) + 1, 

N(u) - N{ji) + 1. (3.21) 

P{jt) changes in the jth bin only by a factor 

/ = (N(js) + l)(iV(j 4 ) + l)/(N( j3 )N(j,)). (3.22) 

Here N(j 3 ) and A^(j 4 ) are the new occupation numbers (after incrementing). R(j') changes 
for each jf > j 

R(j/)^R(j/)+P(j)(f-l) for jf>j. (3.23) 

In addition 



Pint 
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(3.24) 



After making these replacements, one adjusts P(j): 

P(j) - P(j)f. (3.25) 

One then repeats the process as many times as possible choosing new random values of w. 
After very many events the distribution is amplified at two particular values of and p^. 

In the simulation two problems arise. The first is that the enhanced probability to put 
subsequent events in some particular j'th bin is a very mild enhancement at first. Only 
after some huge number of events does the structure lock in on a particular j value. Since 
the number of available fermions in a compact star is of order 10 56 , this is not a problem 
in principle but it does pose technical problems with computers of currently available speed 
and memory. To accelerate the buildup of jet structure, instead of incrementing the selectron 
occupation numbers by one at each throw of the dice, we increment by two. If one increments 
by more than one at each throw, the / of l3.22l has the obvious redefinition. The distributions 
in selectron momentum, polar angle cosine, and azimuthal angle after 600, 000 throws (or 
1.2 million events) is shown in table I. The jagged polar angular distribution shown in figure 
El is an interesting feature of the calculation. 

In a simplified run, (not shown here) where we look at the conversion of two electrons at 
the Fermi surface, although there are no observable jets in the first million events, a clear jet 
structure emerges at 9 • 10 7 events even when selectron occupation numbers are incremented 
by one at each event as is physically required. In this run it is not required that the electron 
pair has zero total momentum in the rest frame of the star. 
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Table 1: Selectron momentum and angular distributions showing the effect of boson en- 
hancement. In this run the occupation numbers are incremented by two at each throw of 
the dice. 
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Figure 3: angular distribution after 600, 000 throws, showing effect of Bose enhancement. In 
this run, the bins are incremented twice at each throw 



The second problem is the following. We have treated the star as having an inexhaustible 
supply of electrons of each momentum in the fermi sea. In practice this number is very large 
but not infinite. Only after all the electrons have been exhausted will the selectron momenta 
add to zero in the rest frame of the star. For the present one can artificially circumvent 
this problem by considering electron to selectron conversion among those electron pairs 
whose center of momentum frame is that of the star. Then we can choose one selectron 
according to the distribution of l3.10l with the other selectron necessarily having the balancing 
momentum. In this case we can rapidly see the jet structure develop incrementing each 
selectron occupation number by one at each throw of the dice. The corresponding selectron 
distribution is shown in figure 0] after only 50, 000 events. 

4 Conclusions 

The results after 600, 000 throws of the dice are shown in table 1. Since the integrated squared 
matrix elements vanish at threshold and the phase space element favors high momenta, the 
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Figure 4: Polar angle distribution for conversion of those electron pairs whose center of 
momentum is that of the star. One quadrant of the angular space is shown with a balancing 
selectron jet in the opposite direction. 



typical selectron momentum is thus somewhat greater than the estimate in [2j based on 
the average energy in the electron Fermi sea. The selectron momentum distributions given 
in table 1 will more or less directly carry over into photon distributions since bosons will 
necessarily fall into the ground state via photon emission. Thus the firm prediction of MeV 
level photons in the burst is a quantitative, and perhaps the primary, result of the current 
paper. A secondary result is that, due to the increased production probability for higher 
photon energy (i.e. selectron kinetic energy in table I), the initial spikes in a gamma ray 
burst would be expected to be of higher average photon energy than later components. This 
is consistent with observations [TT] 

The minimum duration of the burst was fixed in |5j by assuming that the SUSY bubble, 
once nucleated, would grow at the speed of light and that photons in the dense star would 
travel with effective index of refraction one. Both of these assumptions might need to be 
more closely examined. 

In particular, since one should probably regard the bubble surface as a mechanical mem- 
brane, a more reasonable estimate of the burst duration might be given by assuming that the 



bubble expands at the speed of sound rather than the speed of light. Then, if one takes into 
account the range of masses and radii among white dwarfs and if one estimates the bubble 
expansion speed by the speed of sound in a star of corresponding average density, the range 
of predicted burst durations extends well into that of the observations. A detailed study of 
the duration distribution taking into account density inhomogeneities and other effects is a 
high priority subject of future study. 

Finally, we should discuss our expectation that the susy bubble is confined to the dense 
star and does not escape to take over the universe. If one is in the false vacuum of broken 
susy, bubbles of true vacuum (exact susy) are constantly being produced with a steeply 
falling distribution in bubble radius, r. At creation, or at any later stage in its development, 
a bubble of radius r will expand or contract depending on which behavior is energetically 
favorable. The condition for expansion depends on the surface tension, S, of the bubble and 
the energy density in the region immediately outside the bubble. Consider, for example, a 
bubble of exact susy in a larger region of broken susy. Outside the bubble the energy density 
will be e + p where e is the vacuum energy and p is the outside ground state matter density 
if any. If the susy bubble were to make a virtual expansion into an infinitesimally larger 
spherical shell, its ground state energy density would be p s . The difference between p and 
p s is the excitation energy density of the electrons in the broken susy phase. Classically, the 
ground state energy after such a virtual expansion minus the previous ground state energy 
is 

AE = y ((r + 5r) 3 - r 3 ) (p s - p - e) + AnS ((r + 5rf - r 2 ) (4.1) 

or 

AE = -Aixr5r (r(e + Ap) - 25) (4.2) 

where we have put 

Ap = p-p s . (4.3) 
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Classically therefore, the system will find it energetically advantageous to expand if r > 
Similarly, the bubble will contract if its radius is below this density dependent critical 



e+Ap" 

value. A more exact instanton calculation jT^j in vacuo (p = p s = 0) replaces 2S by 35* in 
eq 14.21 One would, therefore, expect the critical radius for a susy bubble to be 

Rc = (4.4) 
e + Ap 

In vacuo or ignoring the effect of the Pauli principle, Ap = 0. In a homogeneous region, if a 
bubble is created at greater than the critical radius, it will expand indefinitely. If however, 
the bubble comes to the boundary of a dense region outside of which p and Ap are zero, the 
critical radius jumps discontinuously to its vacuum value, effectively confining the bubble to 
the high density region. 



Given the indirect hints of supersymmetry from dark matter and accelerator experiments, 
given the positive vacuum energy suggested by the acceleration of the universe, and given the 



persistent suggestion of string theory that the true vacuum of the theory is supersymmetric, 
a phase transition from our broken-SUSY phase to the exact SUSY phase is probably in- 
evitable. There are also good reasons to believe the transition to the true vacuum would be 
catalyzed in dense matter 0j. In addition, given the Bose nature of the final state particles, 
collimation due to stimulated emission is to be expected on physical grounds. Only our 
assumption that the probability of transition is sufficient in compact stars to reproduce the 
rate of gamma ray bursts might be considered speculative. 

It is, however, clear that many points remain to be investigated in the current SUSY 
phase transition picture of gamma ray bursts. Some of these are: 

1. calculation of the hadronic component of SUSY conversion. This will become 
even more essential if the bursts originate in Neutron stars and not in white dwarfs. 

2. calculation of the gamma ray spectrum from radiative SUSY conversion and (not 
independent) bremstrahlung from the final state selectron gas. Since the selectrons are 
bosonic, all of their kinetic energy will emerge as photons. 

3. depletion effects in the Fermi sea 

4. fusion of SUSY (or partially SUSY) nuclei. This could lead to temporary inter- 
ruption of the gamma ray bursts. 

5. polarization effects 

Only after many of these effects have been studied can one attempt to predict the an- 
gular width of jets, the energy spectra, the "light curves", and other features of the rapidly 
growing observational data. As shown here and in |2] the SUSY phase transition model 
provides a framework for discussing other details of this phenomenon that is alternative to 
the more traditional astrophysical approaches. These latter approaches based on relativistic 
ejection of large bodies of neutral matter from black hole accretion disks with the subsequent 
conversion of a large amount of kinetic energy into collimated gamma rays on sub-second 
time scales attempt to describe gamma ray bursts within the boundaries of standard model 
physics but have not as yet led to sharp predictions for typical photon energies or total burst 
energies untied to free parameters in the models. In addition the physical basis of the energy 
release (central engine) or the mechanism for production of collimated gamma rays are not 
yet as fully defined in the accretion models as in the current phase transition model. For 
recent papers marking the current status of traditional astrophysical approaches and giving 
references to earlier work along those lines see and |T3] . 

Another important area that needs study is the possible role of a SUSY phase transition 
in supernova collapse. This is potentially of great interest since it has become apparent 
in recent years that the current standard model of supernova explosion (energy deposited 
by neutrinos in a surrounding shell of matter) is not effective in producing the observed 
explosions ^51 OS]- This problem might also raise questions about the efficiency of vv 
annihilation into an e + e~ cloud available for production of a relativistic fireball to explain 
gamma ray bursts as suggested in [T3| . 



Examples of other new physics suggestions for gamma ray bursts can be found in ^7j 
and references contained therein. Again, it seems that these models might not have the 
same success in predicting the salient features of gamma ray bursts as does the SUSY phase 
transition model. 
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